Spin dependent impurity effects in the 2D frustrated magnetism of NiGa 2 S4 
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Impurity effects on the triangular antiferromagnets Nii_ a! M I Ga2S4 (M = Mn, Fe, Co and Zn) 
are studied. The 2D frozen spin-disordered state of NiGa2S4 is stable against the substitution of 
Zn 2+ (S = 0) and Heisenberg Fe 2+ (S = 2) spins, and exhibits a Independent magnetic specific 
heat, scaled by the Weiss temperature. In contrast, the substitutions with Co 2+ (S = 3/2) spin 
with Ising-like anisotropy and Heisenberg Mn 2+ (S = 5/2) spin induce a conventional spin glass 
phase. From these comparisons, it is suggested that integer size of Heisenberg spins is important to 
stabilize the 2D coherent behavior observed in the frozen spin-disordered state. 

PACS numbers: 75.40.Cx, 75.50.Ee 



Geometrically frustrated and low dimensional magnets 
have attracted considerable interest because of the pos- 
sible emergence of the novel collective phenomena upon 
suppression of conventional magnetic orders. One of the 
most prominent examples of such phenomena is the gap 
formation in ID antiferromagnetic (AF) chains [lj. The 
so-called Haldene gap appears only for chains with inte- 
ger spins. In higher dimensions, the possibility of such a 
2S parity dependent ground state is interesting and in- 
deed have been conjectured by theorists Q , but no exper- 
imental realizations of such an effect has been reported. 

In 2D antiferromagnets (AFMs) based on arrays of tri- 
angles, geometrical frustration in addition to reduced di- 
mensionality may favor collective behavior without spin 
order. One of such exotic phenomena is the quadratic 
temperature dependence in the specific heat without long 
range order, anticipated in a system with 2D coherent 
spin excitations. This phenomenon has been observed 
in an increasing number of AFMs such as the Kagome 
AFMs, SrCrg p Gai2-9pOig [H, deuteronium jarosite 
and more recently in the triangular AFM NiGa2S 4 [5j. 
However, no clear explanation has been made for the 
mechanism to stabilize such a coherent behavior in the 
absence of long range order. 

Here we report the spin dependent impurity effects in 
the 2D coherent behavior in NiGa2S4, the first example 
of the spin 5 = 1 2D AFMs with an exact triangular 
lattice [5J. Our results strongly suggest the 25 parity de- 
pendence in the 2D frustrated magnetism. Interestingly, 
this magnet does not form a conventional 3D AF order 
at least down to 0.08 K in spite of the fact that its AF 
interactions have the energy scale of about 80 K. Instead, 
69 ' 71 Ga nuclear quadrupole resonance and muon exper- 
iments have clarified unusual bulk spin freezing across 
T* = 10 K that has a highly extended critical regime 
down to 2 K Q. Below 2 K, neutron measurements 
have revealed quasi-2D noncollinear correlation, whose 
spin-spin inplane correlation length stays finite around 
7 times the lattice spacing, while the interplane correla- 



tion is so weak that it barely reaches the nearest neighbor 
planes. Moreover, inhomogeneous internal field at the Ga 
site was observed, and the nuclear lattice relaxation rate 
shows nearly cubic temperature dependence below 1 K 
that suggests a 2D magnon-like dispersive mode. Fur- 
ther consistent result has been obtained for the magnetic 
specific heat (Cm) that shows the quadratic temperature 
dependence below about 4 K, suggesting an AF spin- 
wave-like mode in 2D despite the absence of a long range 
AF order. 

Thus, a major question is the origin of the 2D spin- 
wave-like coherent behavior. In order to deepen our un- 
derstanding, we have studied magnetic and nonmagnetic 
impurity effects on the thermodynamic properties of 
NiGa2S4. Our experiments using both single crystals and 
polycrystalline samples of the insulating Nii_ a; M a ;Ga2S4 
(M = Mn, Fe, Co and Zn) indicate that the Fe 2+ (S = 2) 
and Zn 2+ (S = 0) substituted systems, which have 
Heisenberg and integer spins, surprisingly retain the In- 
dependent Cm that is scaled with |#w|- ln contrast, this 
is not the case for the Co and Mn substitutions. The con- 
ventional spin glass phase emerges with the substitution 
of impurities with half-odd integer spins, Ising-like Co 2+ 
(S = 3/2) and Heisenberg Mn 2+ {S = 5/2) spins. This 
suggests that integer size of Heisenberg spins is impor- 
tant to stabilize the 2D spin- wave-like coherent behavior 
observed in the unusual frozen spin disordered state. 

The polycrystalline samples and single crystals of 
Nii_ :c M a ;Ga2S4 are synthesized by solid state reaction 
and chemical vapor transport methods, respectively 0]. 
In order to obtain a homogeneous mixture of Ni and M, 
we first ground the same molar amount of Ni and M 
powder, and repeated the dilutions by adding equal mo- 
lar amount of Ni powder each time until obtaining the 
appropriate M concentration. The amount of M in sin- 
gle crystals was determined by energy-dispersive X-ray 
analyses. Powder X-ray diffraction at room temperature 
on polycrystalline samples indicates single phase isostruc- 
tural to NiGa2S4 at least up to x — 0.1 for Mn, x — 1 for 
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Fe, x = 0.5 for Co and x = 0.3 for Zn. Previous studies 
have clarified that both Zn and Fe atoms occupy the Ni 
site d [§]. For M = Co and Mn, single crystal X-ray 
diffraction experiments were performed, and confirmed 
that Co atoms are indeed located at the Ni site. For Mn, 
the substitution cannot be directly determined because of 
the low concentration. However, the refined cell obtained 
for Mn is isostructural with the trigonal NiGa2S4 and 
consistent with lattice parameters obtained from powder 
diffraction. With Pauling's electrostatic valency principle 
[Io| , one expects Mn atoms to occupy Ni rather than Ga 
site because of the difference in ionic radii (0.66 A and 
0.47 A(Coordination Number, CN = 4) for Mn 2+ and 
Ga 3+ , respectively). Furthermore, the oxidation state as 
obtained from X-ray photoemission spectroscopy (XPS) 
and susceptibility measurements is consistent with the 
Mn 2+ state. DC magnetization M was measured using 
a SQUID magnetometer. The specific heat Cp was mea- 
sured by thermal relaxation method. The lattice con- 
tribution Cl is estimated using the same method as in 
Ref. 0]. Among the impurity ions, Mn and Co have nu- 
clear spins of 5/2 and 7/2, respectively. The Co- nuclear 
quadrupole contribution, C nq = 6.11 x 10 _4 /T 2 (J/mol 
K), was estimated by a standard analysis of the field de- 
pendence. The same procedure was not applicable for 
the Mn case because of the lack of the data points at low 
enough temperatures. Instead, we estimated the nuclear 
quadrupole contribution C nq = 8.55 x 10~ 4 /T 2 (J/mol 
K) by assuming that Mn is located under the same crystal 
field as Co. Finally, Cm was estimated using the relation, 
Cp = Cl + C nq + Cm- 

The high spin (HS) state of Ni 2+ in pure NiGa 2 S4 and 
Mn 2+ has the configuration of (t2 g ) 6 (e g ) 2 and (t2 g ) 3 {e g ) 2 , 
respectively. Both have no orbital degree of freedom, and 
thus their spins should be of Heisenberg type. Fe 2+ has 
the (i2g) 4 (e s ) 2 configuration with HS S — 2, and thus is 
Jahn- Teller (JT) active. A weak rhombohedral distortion 
should split the t2 g orbitals and stabilize one a\ g level 
against the two e' g levels by a JT gap (Ajt), which should 
be of the order of 300 K [9( . Thus, the ground state has 
no orbital degree of freedom, and the spins should be of 
Heisenberg type. For Co 2+ ((t2 g ) 5 (e g ) 2 ), there is residual 
degeneracy of the e' g orbitals, which should lead to Ising- 
like anisotropy of spins, as observed in CoCl2-2H20 [111 ]. 

Figures 1(a) and (b) show the temperature depen- 
dence of the susceptibility \ = M/H and Cp/T for var- 
ious samples. All compounds exhibit spin freezing under 
HqH = 0.01 T as evident from the bifurcation between 
field-cooled (FC) and zero field-cooled (ZFC) sequences. 
Pure NiGa2S4 also exhibits a small bifurcation (~ 6 %) 
below Tf = 8.5 K. In Fig. 1(b), despite the spin freezing, 
Fe 20 % and Zn 10 % substituted compounds surprisingly 
retain the same T 2 -dependence of Cm as NiGa2S4, while 
Cm/T becomes almost constant at low temperatures for 
Mn 10 % and Co 30 % substitution. 

The inset, Fig. 1(c), presents impurity concentration 
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FIG. 1: (color online) Temperature dependence of (a) the 
susceptibility under fioH = 0.01 T, and (b) Cp/T under zero 
field observed by using poly crystalline samples. Impurity con- 
centration dependence of (c) Ax/x at T = 2 K and (d,e,f,g) 
Tf (blue solid square), T p (green solid triangle) and #w (red 
open circle) for all the substituted systems. 



(x) dependence of the ratio (Ax/x) between the size of 
the bifurcation (A%) and the FC susceptibility (%) at 
T = 2 K. All the impurities enhance A%/x with sub- 
stitution. Especially for the Fe substitution, it has a 
maximum value around x = 0.5, suggesting that Ax is 
enhanced by randomness caused by substitution. Thus, 
the bifurcation observed in NiGa2S4 below Tf = 8.5 K 
may be attributed to freezing of minority spins following 
the bulk spin freezing at T* — 10 K. While the increase 
in Ax with substitution is similar to what has been ob- 
served in conventional spin glass systems, a clear differ- 
ence appears in the change in Tf. Figures 1(d), (e), (f) 
and (g) show ^-dependence of Tf , peak temperature T p of 
Cm/T, and &w, where Tf is determined to be the bifur- 
cation point under 0.01 T and 6\y is by the Curie- Weiss 
(CW) analysis using the formula x — C/(T — W ). Nor- 
mally, for conventional spin glasses, Tf increases with the 
impurity concentration. However, Tf for all the impurity 
substitutions except Fe decreases in the small concentra- 
tion regime (x < 0.1). Moreover, for the Fe substituted 
system, Tf decreases with Ni substitution for FeGa2S4, 
and this is not consistent with conventional spin glass 
behavior. Instead, the Fe, Zn and Mn substituted sys- 
tems show a scaling behavior, |0vv| ~ 10Tf, indicating 
that Tf is determined by #w- However, this is not the 
case only for the Co substituted systems. Meanwhile, T p 
roughly scales with Tf but always slightly above Tf , sug- 
gesting the broad peak in Cm/T is associated with the 
spin freezing. Only one exception for this case is the Fe 
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FIG. 2: (color online) (a,b) Effective moments and spin size 
of M in the polycrystalline samples, (c) Temperature depen- 
dence of the ratio between \ab and Xc- 



substituted system. This is because of the existence of 
Ajt, which will be discussed later. 

Effective moments obtained by the CW fitting to the 
susceptibility data of polycrystalline samples are plotted 
in Fig. 2(a) and (b). The results are consistent with the 
value (solid line), 2y2(l — x) + 2-y/5i mp (5j mp + l)x, es- 
timated for HS state, where 5; mp is the size of impurity 
spin. Furthermore, a systematic change of the lattice 
constants with substitution [?J and recent XPS measure- 
ments by using single crystals [13] also suggest that all 
the M ions have 2+ valence and are in the HS states. In 
Fig. 2 (a), an additional enhancement of p c g was found in 
the Co substituted system, which should come from the 
orbital contribution due to the degeneracy of e' g orbitals. 

Figure 2(c) shows anisotropy determined by the ratio 
between the inplane and interplane susceptibility (xab 
and Xc) under 0.01 T FC sequence. NiGa2S4, FeGa 2 S4, 
and the Mn, Zn substituted systems have 1 < Xab/Xc < 
1.4, which indicates that spins are of Heisenberg type 
with weak XY anisotropy, while the Co substituted sys- 
tem shows Ising-like anisotropy with Xab/Xc down to 0.5. 
The observed anisotropy is consistent with the single-ion 
anisotropy expected for each impurity ion. 

For states that have a single energy scale of |#\y|, the 
molar entropy (5m) should take the form 5m = i?ln(25+ 
l)f(T/\(hjv\) including its spin size (5). By taking the 
derivative of this, one obtains the following relation, 
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In order to check the normalization with |#w|> 
C M \0w\/(TR\n(2S + 1)) vs. T/|0 W | is plotted in Fig.3, 
where the total molar spin entropy in the vertical axis is 
given by R{(1 — x) In 3 + xln(25i mp + 1)}. Figures 3(a) 
and (b) are for the Zn and Fe substituted systems. The 
Zn substituted system exhibits the T 2 -dependent behav- 
ior of Cm up to x = 0.3. The substitution leads to the 
abrupt depression of the peak in Cm/T by the substitu- 
tion of less than 1 %, indicating our obtained sample of 
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FIG. 3: (color online) T/\9w\ dependence of the normalized 
Cm for (a) Zn, (b) Fe, (c) Co and (d) Mn substitutions. In- 
dependence of Cm/T for (e) Co and (f) Mn substitutions. 
Solid lines indicate its linear extrapolation below 0.6 K. 



NiGa2S4 is near the clean limit. The initial slope remain 
constant, indicating that there exists a normalized func- 
tion f(T/\6w\) and the T 2 -coefficient of Cm is propor- 
tional to |6>w|~ 2 - This (T/|# w |) 2 form of the specific heat 
suggests an existence of a 2D spin-wave-like mode whose 
spin stiffness is proportional to #w- The two dimension- 
ality is also consistent with the fact that 6\y should be 
determined not by the negligibly small interplane cou- 
pling, but by the inplane exchange interactions. For the 
Fe substitution, all the curves have the well-scaled initial 
slope up to T/\8-w\ ~ 0.05, above which, however, they 
start deviating from each other. This may be due to the 
effect of the JT gap A JT of Fe 2 + ions. Given Ajt ~ 300 
K [9(, the same order of magnitude as |#w|, the Schottky 
tail due to Ajt might affect the T-dependence of Cm/T 
and its peak formation at ~ T p . 

In Fig. 3(c), the normalized Cm for the Co substitu- 
tion is shown. It is clear that no curves overlap each 
other, indicating that Co spins with Ising anisotropy 
completely suppress the scaling behavior of Cm/T. More- 
over, 7 which is estimated by the linear extrapolation of 
Cm/T curve below 0.6 K (Fig. 3(e)), starts appearing at 
x > 0.1 and increases with the Co content. The anomaly 
at x — 0.1 is also seen in Tf, T p (Fig. 1(f)) and Ax/x at 2 
K (Fig. 1(c)), pointing toward the existence of a crossover 
at x ~ 0.1. At x < 0.1, the system exhibits the behav- 
iors not consistent with the conventional spin glass, i.e. 
Tf decreases against impurities, and 7 keeps its value 
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FIG. 4: (color online) T/|0 W | dependence of C M /T 2 - A for 
the (a) Zn, (b) Fe, (c) Co and (d) Mn substituted systems. A 
is determined by using Eq.(2). 



~ with Cm ~ T 2 -behavior. At x > 0.1, on the other 
hand, the system shows conventional spin glass behavior 
in the sense that it shows finite 7 and that both 7 and 
Tf increase with the substitution. These results indicate 
that the ground state of NiGa2S4 is different from that 
of conventional spin glass state. 

Finally, for Mn substitution, even with the relatively 
low amount (x < 0.1) of impurity, no scaling is found 
in the normalized Cm/T (Fig. 3(d)), but a systematic 
change from the Independence of Cm (Fig. 3(f)). In order 
to show this explicitly, we plot in Fig. 4 the deviation in 
Cm from the Independence estimated using the scaling 
law of Eq.(l) as follows. Suppose Cm = AT 2 behav- 
ior is well scaled by a single energy scale of |#w|j A for 
Nii_ a; M I ;Ga2S4 should take the form, 

. E{(l-x)m3 + xm(2S imp + l)} 

A = a U{ — P5 ! W 



\8w\ 2 



leading /'(T/|0 W |) = aoT/\6 w \ in Eq.(l). Here, the co- 
efficient ao = 19.4 is uniquely determined by the T 2 - 
coefficient of Cm for pure NiGa2S4. |0w| is determined 
by the CW analysis of x(T), thus there is no free param- 
eter in this analysis. The vertical axis in Fig. 4 presents 
Cm/T 2 — A, using A estimated by Eq.(2). If the scaling 
works, the vertical axis in Fig. 4 should take constant zero 
value. Indeed, this is the case for the Zn and Fe substi- 
tution (Figs. 4(a) and (b)). However, the Mn and Co 
substituted systems do not hold the Cm ~ T 2 behavior 
(Figs. 4(c) and (d)), but have a finite 7. These indicate 
that a conventional spin glass phase is stabilized by the 
Mn substitution, as for the Co case with x(Co) > 0.1. 

To summarize the experimental results, a conventional 
spin glass phase emerges with the substitution of mag- 
netic impurities with half-odd integer spins: Ising-like Co 
spins and Heisenberg Mn spins. Such conventional spin 



glass behavior is expected, given that the randomness 
due to impurities and geometrical frustration inherent in 
the triangular lattice. Surprisingly, however, the Fe and 
Zn substituted systems which have Heisenberg and in- 
teger spins, retain the Independent Cm, that is scaled 
with |#w|- Moreover, Tf also scales with 6w These sug- 
gests that the unusual spin frozen state of NiGa2S4 is 
not conventional spin glass, and its 2D coherent behavior 
represented by Cm ~ (T/|#\y|) 2 appears only for integer 
Heisenberg spins. 

Two types of theoretical proposals have been made for 
the origin of the 2D coherent behavior in NiGa2S4. First 
one is the spin nematic state that describes unconven- 
tional long-range order by magnetic quadrupoles without 
long-range two-spin correlation. In this state, the site av- 
erage of spin is zero [3, 0, EH, EH , which is inconsistent 
with observed results such as the internal field. Second 
one is a Kosterlitz-Thouless (KT) type phase driven by 
Z2 vortices [l7j |. In this case, the anomaly in Cm is due 
to the condensation of Z2 vortices that are formed by 
vector spin chirality. It is likely that Ising type impu- 
rity spins destabilize such crossover as in our experiment 
because vector spin chirality can be only defined for the 
Heisenberg spins. A recent theory [181 ] based on quantum 
calculations that treats the low-T region below the KT 
type crossover predicts the (T/|6\y Independent specific 
heat and the constant susceptibility as T — > K, which 
is consistent with our observations. 

However, this Z2 vortex mechanism alone cannot ex- 
plain the 2S parity dependence observed in the present 
experiments. The distinct behavior between integer and 
half-integer spins may be understood as a quantum effect, 
for example, through quadratic form of spin exc hang e in- 
teractions, as known in Haldane systems in ID 19]. In- 
deed such a quadratic form of spin interactions has been 
discussed to explain the magnetism of the related com- 
pound MS2 [201 ]. This type of quantum effect might be 
also the origin of the finite inplane correlation length. 
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